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ABSTRACT 

Data on the fluctuations in depth of maximum development of 
cosmic ray air showers, corrected for the effects of mixed 
primary composition and shower development fluctuations, 
yield values of the inelastic proton-air cross section for la- 
boratory energies in the range 10 8 -10 10 GeV. From these val- 
ues of 0par corresponding values of the proton-proton total 
cross section are derived by means of Glauber theory and geo- 
metrical scaling. The resulting values of Opp are inconsis- 
tent with a well known ln 2 s extrapolation of ISR data which is 
consistent with SPS data; they indicate a less rapid rate of 
increase in the interval 540 < /s < 10 5 GeV. 

1. Introduction. The total inelastic proton-air cross section a pa , which 
is interesting in relation to the theory of hadron— hadron and hadron— nuc- 
leus interactions, can be derived from data on cosmic ray air showers. 
Cosmic ray evidence of an increase with increasing energy in c?pp, the to- 
tal proton-proton cross section, has been confirmed, first at ISR ener- 
gies and more recently with the CERN SPS collider, 1 " 4 in the interval 
10^< /s < 540 GeV the increase, amounting to some 75%, agrees well with a 
In s dependence but is compatible with a slower rate. 3 A recent cosmic 
ray result for 2*10 < /s < 2*10 4 GeV favors ln 2 s, 9 while another result 

of this kind, for /s = 3»10 4 GeV, favors a slower increase. 6 In case of 
the cosmic ray results, CFpp is derived from Cp a using Glauber theory. 

M Y purpose here is to draw attention to another body of air shower 
evidence extending to even higher energies (10 4 < /s < 1.6- 10 5 GeV) which 
also bears on this question. 19 The evidence consists in part of data on 
the fluctuation of shower elongation x m , and in part of data on the pri- 
mary composition. These are combined using an expression for {x^J as a 
function of <lnA>, {lnA} and {x^p}, where A is the primary mass number, 
<lnA>, {lnA} are the mean value and dispersion about the mean of lnA for 
fix^d primary energy, {x m ^p} is the dispersion of x m for proton-initiated 
showers of a fixed energy (allowing for fluctuations in shower origin and 
development) , and {x m } is the observed dispersion of x m , corrected for 
reception fluctuations (instrumental dispersion). 7 Using experimental 
evidence on <lnA> and {lnA} from other air shower observations, the cor- 
rection for contamination of the primary beam with heavy nuclei is shown 
to be small; { x m,p^ = (x^. Cascade simulations have shown that {x mf p} 
is proportional to Ap a , the proton-air mean free path, with a proportion- 
ality constant equal to 1.4. From Ap a one obtains 0p a ; Jpp follows by 
application of Glauber theory. 

2* P^iraary Composition Effect. I have shown elsewhere that if y is an 
observable such as x m , whose mean value for proton showers can be graphed 
vs InE as a straight line (over reasonably wide intervals, say 2 decades 
or more) , and if the line width (fluctuation of y) for showers initiated 
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by various nuclei is approximated by {y(A)} = {y(l)}(l - klnA) , where k 
is constant, then 

{ y } 2 = {y (1) } 2 [1 - 2k<lnA> + k 2 (<lnA> 2 + {InA} 2 ) ] + b 2 {lnA} 2 , (1) 

where b is the slope of the straight-line graph. 7 This follows from ele- 
mentary statistical theory and the superposition principle of Peters. 

For y = x mf b is just D S/ p, the proton shower elongation rate. With 
these substitutions I solve (1) for {x m> p}, obtaining 


r 1 - (D { InA}/ {x }) 

K } = {x } — 2 — 

m,p m l 1 - 2k<lnA> + k 2 (<lnA> + {InA} ) 


The right hand side of (2) contains, besides the observed quantity re- 
quiring correction, 2 parameters which describe phenomenologically the 
relevant features of air showers (D Gf p and k) , and 2 parameters which de- 
scribe the composition of the incoming cosmic ray beam (<lnA> and {InA}) • 

The fluctuation of x m has been studied experimentally for more than 
10 years using the giant air shower r.rray at Haverah Park, 8 ' 9 and almost 
as long using a similar size array at Yakutsk. 10 It has been measured 
independently at Dugway using an array of atmospheric Cerenkov light de- 
tectors f 11 and at Akeno using data on the electron lateral distribution. 
Results are given in Table 1. 

I have shown elsewhere that under a wide range of assumptions about 
hadron-nucleus interactions D e , p cannot be appreciably greater than t Q , 
the radiation length. This insures that the numerator in (2) is approx- 
imately 1, so the exact value of D e/ p is unimportant. For energy inde- 
pendent composition, D G/P = D e (observed elongation rate). I adopt D e ,p 
= 28 g/cm 2 , practically * the experimental value of D e at these energies. 4 
As for the parameter k, Monte Carlo cascade simulations show that its 
value lies in the range 0.15 ± .05, depending on the choice of model. 

The exact value is unimportant if <lnA> is less than 1. 

Data on <lnA> and {InA} at very high energies are summarized in Ref. 
7 (see also conference paper OG5.4-4),. They can be combined, as in Fig. 
1, by plotting {InA} vs <lnA> in rectangular 
coordinates. Note that for any pure compo- 
sition {InA} =0. It is assumed that the 
possibilities range from <lnA> = 0 (pure 
protons) to <lnA> = 4 (pure Fe) . The maxi- 
mum dispersion for a given mean occurs for a 
mixture of protons and Fe nuclei. When the 
scales are chosen as I have done the locus of 
binary proton— Fe mixtures is a semicircle. 

The lower boundary of the diagram consists of 
smaller semicircles corresponding to other 
binary mixtures. The region inside the 
boundary represents all possible composi- 
tions. The observed equal-energy composi- 
tion in the low energy region (10 2 -10 5 GeV) Fig, 1. Cosmic ray com - 

is indicated by the heavy dashed line. The position diagram: 2nd vs 

compositions experimentally allowed at very 1st moment of the InA dis - 
high energies (E > 10 8 GeV) are those within tribution, explained in 
the shaded region. The heavy dotted line the accompanying text. 
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Table 1 . Summary of {x^} data 


represents a detailed model given in 
Ref, 7, which describes in a plausible 
manner the high energy region and the 
transition from lower energies. The 
correction factors {x mf p}/{x m } for this 
model to be used with Table 1 according 
to (2) range in value from 0.91 to 0.95; 
thus they differ only slightly from 1 
(no correction) . These corrections have 
been applied in obtaining the present 
results, shown in Figures 2 and 3. 

Also shown in Fig. 1 are contour 
lines corresponding to a typical value 
of {x^} (= 66 g/cm 2 ) and the indicated 
values of the correction factor. By ex- 
amining these contours one sees that 
correction factors for any other allowed 
model will differ only slightly from 
those used here. 

3. Development Fluctuations. The next 
step is to convert from {x m ,p} to a pa 
allowing for shower development fluctua- 
tions. These arise primarily because of 
the leading particle effect. Thus one 
expects them to scale as Ap a , so there 
will be a proportionality: {x m ^p} = KA pa 
where K is a constant. On the basis of Monte Carlo cascade simulations 
by Walker and Watson^ I adopt the value 1.40 for this constant. The re- 
sulting values for a pa are given in Figure 2. 

There is of course some uncertainty in the value of K, due mainly to 
uncertainty about the elasticity of very high energy hadron-air colli- 
sions, denoted here by ri. The result K = 1.4 rests on a conventional as- 
sumption that n is distributed uniformly between 0 and 1 with a mean val- 
ue of 0.5. It has been pointed out by L.W. 

Jones that according to accelerator data 
the value of <ri> for proton collisions with 
light nuclei is appreciably less, about 
0.3. 16 The dependence of K on <T)> can be 
estimated by assuming 1) that x av of the 
partial electronic cascade originating 
from each collision of the leading nucleus 
is perfectly correlated with the depth of 
that collision, and 2) that {x m ,p} = 

^ x av,p) (neglecting fluctuations in n) . 

It is then easily shown that 

(3) 
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total of 464 events 


00 _ 1 -i. 

K = j (1-n) n J r - 


(i-n) 


For n = 0.5 one obtains K s 1.41; for n = 
0.3 the result is K = 1.20. Thus the ef- 
fect of assuming a smaller value of <n> is 
to reduce the apparent cross sections o r 
and Opp at air shower energies. 
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the symbols are the cor- 
responding references . 




4, Results. Examining Fig. 2 one notes good 

agreement between the present results and one 
published recently by the Fly's Eye group. 6 
The Fly's Eye result is derived from the de- 
crement of the x m distribution, rather than 
the dispersion, using a proportionality be- 
tween and Ap a . The value of the constant 

is taken as 1.6 ± 10% from simulations 1 5 si- 
milar to those of Walker and Watson. The 
present results disagree with those of Ref. 

5, based on the zenith angle distribution of 
generally smaller air showers for fixed N e 
and Ny. Here also the relation between a 
measured attenuation length and Ap a is found 
by means of simulations. There is an approx- 
imate proportionality, the value of the con- 
stant being 1.45-1.55. The value of <n> for 
the Akeno model is 0.44. A smaller value 
would result in a smaller value of the pro- 
portionality constant and a smaller apparent 
rise in cross sections. However it would not 
remove the discrepancy since a change in <n> 
will have a similar effect on results from 

A xm and **m>- 

Fig. 3 shows the result of converting 
the new a pa values to Gpp. This has been 
done the same way as in Ref. 6, using Glauber 
theory with an assumption that the nuclear 
slope parameter has the same energy dependence 
as Gpp (geometrical scaling) . In parametric form, the relation assumed 
is cr pa =5 24.0 Gpp° * 6 48 (cross sections in mb). Shown for comparison are 
the SPS results, 3 ' 4 and well known extrapolations by Block and Cahn, 
based on ISR data, which successfully predicted those results. 17 Also 
shown is a recent prediction by Bourrely et al., for the so-called impact 
picture. 18 The lower Block-Cahn curve is the best fit with a ^ 0, where 
parameter a takes on small positive values to allow for deviations from 
Froissart bound form. The upper curve is for a = 0. 
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Fig. 3 . Results on o pp 
(same symbols as Fig . 2 ) . 
The lower solid curve , 
from Ref . 17, is a best 
fit to ISR data allowing 
deviation from the Frois- 
sart bound. The upper 
solid curve (same Ref.) 
is a best fit with Frois- 
sart bound form . The 
dashed curve is for the 
9 impact picture 1 . 18 




